Black Holes at the IceCube Neutrino Telescope 
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If the fundamental Planck scale is about a TeV and the cosmic neutrino flux is at the Waxman- 
Bahcall level, quantum black holes are created daily in the Antarctic ice-cap. We re-examine the 
prospects for observing such black holes with the IceCube neutrino-detection experiment. To this 
end, we first revise the black hole production rate by incorporating the effects of inelasticty, i.e., 
the energy radiated in gravitational waves by the multipole moments of the incoming shock waves. 
After that we study in detail the process of Hawking evaporation accounting for the black hole's 
large momentum in the lab system. We derive the energy spectrum of the Planckian cloud which 
is swept forward with a large, O(10 6 ), Lorentz factor. (It is noteworthy that the boosted thermal 
spectrum is also relevant for the study of near-extremal supersymmetric black holes, which could 
be copiously produced at the LHC.) In the semiclassical regime, we estimate the average energy 
of the boosted particles to be less than 20% the energy of the ^-progenitor. Armed with such a 
constraint, we determine the discovery reach of IceCube by tagging on "soft" (relative to what one 
would expect from charged current standard model processes) muons escaping the electromagnetic 
shower bubble produced by the black hole's light descendants. The statistically significant 5<r excess 
extends up to a quantum gravity scale ~ 1.3 TeV. 
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I. GENERAL IDEA 

Over the past few years it has become evident that 
a promising approach towards reconciling the apparent 
mismatch of the fundamental scales of particle physics 
and gravity is to modify the short distance behavior of 
gravity at scales much larger than the Planck length, 
Zpi ~ 10~ 35 m. The key premise of such an approach 
entails that the weakness of gravity is indeed evidence of 
extra spatial compactified dimensions [IJ. This is possi- 
ble because standard model (SM) fields are confined to a 
4-dimensional world (corresponding to our apparent uni- 
verse) and only gravity spills into the higher dimensional 
spacetime bulk, without conflicting with experimental 
bounds 0]. Therefore, if this new approach is correct, 
gravity is not intrinsically weak, but of course appears 
weak at relatively large distances of common experience 
because its effects are diluted by propagation in the extra 
dimensions. The distance at which the gravitational and 
electromagnetic forces might have equal strength is un- 
known, but a particularly interesting possibility is that 
it is roughly at 10 -19 m, the distance at which electro- 
magnetic and weak forces are known to unify to form the 
electroweak force. This would imply a fundamental D- 
dimensional Planck mass, Mr> ~ M\y ~ 1 TeV, consid- 
erably smaller than the macroscopic 4-dimensional value, 
Mpi - 10 19 GeV. 

If nature gracefully picked a sufficiently low-scale grav- 
ity, the first evidence for it would likely be the observa- 
tion of microscopic black holes (BHs) produced in parti- 
cle collisions Q. According to Thome's hoop conjecture, 
a BH forms in a two-particle collision when and only 
when the impact parameter is smaller than the radius of 
a Schwarzschild BH of mass equal to the total center-of- 
mass energy Q . Subsequent to formation a TeV-scale BH 
will promptly decay via thermal Hawking radiation [l7j 
(for M D = 1 TeV, the lifetime of a BH of mass 10 TeV is 



less than 10~ 25 s) into observable quanta [5[. Although 
the BH production cross section, 0(M^}), is about 5 
orders of magnitude smaller than QCD cross sections, 
0(Aqq D ), in two well-known papers [f| 0] it was pro- 
posed that such BHs could be produced copiously at the 
LHC, and that these spectacular events could be easily 
filtered out of the QCD background. This is possible by 
triggering on BH events with prompt charged leptons and 
photons, each carrying hundreds of GeV of energy. 

Cosmic ray collisions, with center-of-mass energies 
ranging up to 10 5 GeV, certainly produce BHs if the 
LHC does. The question is, can they be detected? Most 
cosmic rays are protons, which generally collide with 
hadrons in the upper atmosphere, producing cascading 
showers which eventually reach the Earth's surface. At 
energies of interest, however, the cosmic ray luminos- 
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cm 2 s , is about 50 orders of mag- 



nitude smaller than the LHC luminosity, thus making 
it futile to hunt for BHs in baryonic cosmic rays. On 
the other hand, neutrino interaction lengths are still far 
larger than the Earth's atmospheric depth, although they 
would be greatly reduced by the cross section for BH 
production Q. Cosmic neutrinos therefore would pro- 
duce BHs with roughly equal probability at any point 
in the atmosphere. As a result, the light descendants of 
the BH may initiate low-altitude, quasi-horizontal show- 
ers at rates significantly higher than SM predictions 0]. 
Because of these considerations the atmosphere provides 
a buffer against contamination by mismeasured baryons, 
for which the electromagnetic channel is filtered out. 

Neutrinos that traverse the atmosphere unscathed can 
produce BHs through interactions in the Antarctic ice- 
cap and be detected by the IceCube neutrino tele- 
scope [HI- This telescope, which is currently being de- 
ployed near the Amundsen-Scott station, comprises a 
cubic-kilometer of ultra-clear ice about a mile below the 
South Pole surface, instrumented with long strings of sen- 
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sitive photon detectors which record light produced when 
neutrinos interact in the ice [111 ] . The In- ice array is com- 
plemented by IceTop, a surface air shower detector con- 
sisting of a set of 160 frozen water tanks, which serves as 
a veto for atmospheric muon background. At the same 
time, the energy deposited by tagged muon bundles in air 
shower cores becomes an external source of energy cal- 
ibration. Altogether, the expected energy resolution of 
the experiment is about ±0.1 on a logio scale. Moreover, 
the energy reconstruction is optimized for neutrino en- 
ergy E v > 10 6 GeV, allowing sufficient precision (±0.2 on 
a logio scale) to separately assign the energy fraction for 
emergent muons in neutrino interactions. Because of this, 
the inelasticity distribution of events becomes a unique 
tool for SM background rejection, providing powerful dis- 
crimination of resonant processes (T3 |. In this work we 
re-examine the prospects for discovering BH resonances 
at IceCube, by tagging on "soft" (relative to what one 
would expect from charged current SM processes) muons 
escaping the electromagnetic shower bubble triggered by 
the BH explosion. 

The paper is organized as follows. In Sec.[ll]we update 
the semi-classical BH production cross section consider- 
ing a new estimate [l3[ of the energy radiated in gravita- 
tional waves by the multipole moments of the incoming 
shock waves. This is followed by a detailed discussion 
of Hawking evaporation taking into account that cosmic 
neutrinos produce BHs with large momenta in the lab 
system. Specifically, we derive the energy distribution 
of the Planckian cloud which is boosted in the forward 
direction with a large Lorentz factor. Armed with this 
distribution, we estimate the average energy of the BH 
light descendants to be less than 20% the energy of the 
z^-progenitor. The Hawking radiated muons then pro- 
vide a very clean signal with negligible SM background, 
as the production of "soft" leptons in charged current 
(CC) interactions occurs at a much smaller rate than BH 
production. Our results for event rates and discovery 
reach are presented in Sec. IIII1 Conclusions are given in 
Sec. lIVl In this last section we also entertain the possibil- 
ity of producing TeV-scale near-extremal BHs in particle 
collisions. 



II. BH PRODUCTION AND EVAPORATION 

Analytic and numerical studies have revealed that 
gravitational collapse takes place at sufficiently high en- 
ergies and small impact parameters, as conjectured years 
ago by Thorne [J] . In the case of 4-dimensional head-on 
collision s ww . as well as those with non-zero impact pa- 
rameter |l5| , a horizon forms when and only when a mass 
is compacted into a hoop whose circumference in every 
direction is less than 2ir times its Schwarzschild radius up 
to a factor of order 1. In the D-dimensional scenario the 
Schwarzschild radius still characterizes the maximum im- 
pact parameter for horizon formation [16]. In the course 
of collapse, a certain amount of energy is radiated in 



gravitational waves by the multipole moments of the in- 
coming shock waves [1J], leaving a fraction y = Mbh/Vs 
available for Hawking evaporation [llj. Here, Mbh is a 
lower bound on the final mass of the BH and y/J is the 
center-of-mass energy of the colliding particles, taken as 
partons. This ratio depends on the impact parameter of 
the collision, as well as on the dimensionality of space- 
time [3. 

Of course, this work is purely in the framework of 
classical general relativity, which is valid only for suf- 
ficiently massive BHs, A/bh ^> Mp. For masses close 
to Md, gravity becomes strong and the classical descrip- 
tion can no longer be trusted. Hence, it is important to 
impose a lower cutoff on the mass of microscopic BHs 
for which the simple semiclassical arguments can rea- 
sonably be expected to hold. Following we define 
Xmin = Mg^/Mo = 3, where Mgg 1 is the smallest BH 
mass for which we trust the semiclassical approximation. 

String theory provides a promising route for under- 
standing the regime of strong quantum gravity and in 
particular for computing cross sections at energies close 
to the Planck scale [20] ■ Therefore, the ensuing dis- 
cussion will be framed in the context of string theory. 
To be specific we will consider an embedding of a 10- 
dimensional low-energy scale gravity scenario within the 
context of SO (32) Type I superstring theory, where gauge 
and charged SM fields can be identified with open strings 
localized on a 3-brane and the gravitational sector con- 
sists of closed strings that propagate freely in the internal 
dimensions of the universe [211 ]. After compactification 
on T 6 down to four dimensions, Mpi is related to the 
string scale, M s , and the string coupling constant, g s , by 
Mpj = (2-7T r c ) 6 Mg/<7g , where r c is the compactification 
radius. Within this framework, the problem of avoid- 
ing fast baryon decay [22[ or lepton flavor violation [23| 
is shifted to the examination of symmetries [24| in the 
underlying string theory which would suppress the ap- 
propriate non-renormalizable operators at low energies. 
Nevertheless, it is important to stress that for x m [ n > 3 
the typical decay involves a large number of particles. 
Therefore, though these symmetries constrain the decay 
of the BH, throughout this paper we ignore the con- 
straints imposed by the few conservation laws and we 
assume that BHs decay with roughly equal probability 
into all SM particles. From now on we set D — 10. 

The inclusive production of BHs proceeds through dif- 
ferent final states for different classical impact param- 
eters b [l||. These final states are characterized by 
the fraction y(z) of the initial parton center-of-mass en- 
ergy, vS) which is trapped within the horizon. Here, 
z = b/b max , and & ma x = V~F r s (s/§) is the maximum im- 
pact parameter for collapse, where 
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is the radius of a Schwarzschild BH in 10-dimensions [2f 
and F is a form factor. 
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A bound on the inelasticity and the form factor can 
be obtained by studying the formation of an aparent 
horizon, which (because of cosmic censorship [261] ) gu- 
rantees the formation of a BH event horizon [27| . Such a 
study can be easily accomplished by modeling the incom- 
ing partons as two Aichelburg-Sexl [28| shock waves (i.e, 
by boosting the Schwarzschild solution to the speed of 
light at fixed energy). The scattering of partons is then 
simulated through the superposition of two shock waves 
coming from opposite directions, such that their union 
defines a closed trapped surface which provides a lower 
bound on Mbh and 6 ma x [Uj]- This lower bound, how- 
ever, depends on the slice used to determine the apparent 
horizon, and becomes larger if the apparent horizon is 
taken on the future light cone of the collision plane [HI . 
This is because it is possible that for a given impact pa- 
rameter an apparent horizon is not yet formed on the 
so-called "old slice" , but arises by the time a later "new 
slice" is reached. In our calculations we consider the es- 
timates of y and F in both the "old" [l8| and "new" [l3T ] 
slices. 

The y dependance complicates the parton model cal- 
culation, since the production of a BH of mass Mbh re- 
quires that s be M^- li /y 2 (z), thus requiring the lower 
cutoff on parton momentum fraction to be a function of 
impact parameter [29( . Because of the complexity of the 
final state, we assume that amplitude intereference ef- 
fects can be ignored and we take the vN cross section 
as an impact parameter-weighted average over parton 
cross sections, with the lower parton fractional momen- 
tum cutoff determined by x mm . This gives a lower bound 
X = (x ' m i n M io) 2 /[y 2 (z)s] on the parton momentum frac- 
tion x, where a/s is the center-of-mass energy of the vN 
collision. All in all, the vN — ► BH cross section reads [3(| 

a= f 2zdz [ dxFm 2 s {V~s,M w ) 5~)/i(s,Q) , (2) 

Jo J X \ 

where s = xs = 2xm^E u , i labels parton species, and 
the fi(x,Q) are parton distribution functions (pdfs). 

The choice of the momentum transfer Q is governed 
by considering the time or distance scale probed by the 
interaction. Roughly speaking, the formation of a well- 
defined horizon occurs when the colliding particles are at 
a distance ~ r s apart. This has led to the advocacy of 



the choice Q 



[311 ]. which has the advantage of a 



sensible limit at very high energies. However, the dual 
resonance picture of string theory [32] would suggest a 
choice Q ~ ■</§. Fortunately, as noted elsewhere [33| . 
the BH production cross section is largely insensitive to 
the details of the choice of Q. In what follows we use the 
CTEQ6D pdfs H with Q = minjrj 1 , 10 TeV}. In Fig.Q] 
we show the BH production cross section for x min = 3 
and Mio = 1 TeV. 

BHs produced in particle collisions have non- vanishing 
angular momenta determined by the impact parameter 
of the incoming partons. Since b is only non-zero along 
the brane directions, the angular momentum lies within 
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FIG. 1: BH production cross section for x m i n = 3 and Mio = 
1 TeV. The solid line indicates the result obtained using the 
new estimates of F and y given in , whereas the dot-dashed 
line indicates the result obtained assuming the old values from 
Ref. IH. 



the brane. Moreover, the initial horizon is likely to be 
very asymmetric as only gravity spills into the compact- 
ified dimensions. Therefore, the excited BH state carries 
additional hair corresponding to multipole moments for 
the distribution of gauge charges and energy momentum 
within the asymmetric configuration. The decay of an 
excited spinning BH state follows several stages. The 
initial configuration looses hair associated with multi- 
pole moments in a balding phase through emission of 
gravitational waves. In addition, gauge charges inherited 
from the initial state partons are discharged within this 
phase via Schwinger emission. The subsequent spinning 
BH evaporates in a two-step process: a short spin-down 
phase in which angular momentum is shed [35| . followed 
by a long Schwarzschild phase of semi-classical Hawking 
radiation. 

In the rest frame of the Schwarzschild BH, both the 
average number [FH ] and the probability distribution of 
the number [3(| of outgoing particles in each mode obey 
a thermal spectrum. In 10-dimensions, the emission rate 
per degree of particle freedom i of particles of spin s with 
initial total energy between (u>, ui + du>) is found to be [13] 
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is the instantaneous Hawking temperature, 
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is the volume of a unit (d — 3)-sphere, 
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(6) 



is the instantaneous Schwarzschild radius of mass M, 
and <7 s (w) is the greybody absorption area due to the 
backscattering of part of the outgoing radiation of fre- 
quency u into the BH (a.k.a. the greybody factor) [38j |. 
Recall that SM fields live on a 3-brane (d — 4), while 
gravitons inhabit the entire spacetime (d = 10). The 
prevalent energies of the decay quanta are of 0(T ~ 1/r), 
resulting in s-wave dominance of the final state. Indeed, 
as the total angular momentum number of the emitted 
field increases, a s {u>) rapidly gets suppressed [39j]. In the 
low energy limit, wr « 1, higher-order terms are sup- 
pressed by a factor of 3(wr) -2 for fermions and by a 
factor of 25(w r)~ 2 for gauge bosons. For an average par- 
ticle energy (us) of C(r _1 ), higher partial waves also get 
suppressed, although by a smaller factor. This strongly 
suggests that the BH is sensitive only to the radial coor- 
dinate and does not make use of the extra angular modes 
available in the internal space Actually, a recent de- 
tailed analysis [40| has explicitly shown that the relative 
emission rate of SM particles and the 10-dimensional bulk 
graviton is roughly 92:5. This implies that the power lost 
in the bulk is less than 15% of Mbh, largely favoring the 
dominance of visible decay. Therefore, in what follows 
we assume the Hawking evaporation process to be domi- 
nated by SM brane modes and we neglect graviton emis- 
sion during the Schwarzschild phase. With this in mind, 
the average total emission rate for particle species i is, 



d(N) 
dt 
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where Cj is the number of internal degrees of freedom of 
particle species i, gi = 1 (3/4) for bosons (fermions), 
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The rate of change of the BH mass in the evaporation 
process is 

dM 

~ 2tt 



dt 



r(E C(4)T(4)r 2 T 4 , 



cvap 
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where /» = 1 (7/8) for bosons (fermions) and 
<j s {uj) to 3 dio 
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Dividing Eq. ([7]) by Eq. ^ and integrating, one obtains 
a compact expression for the average multiplicity 

(N) = \ P [8 ^r(9/2)] V7 [^l 8/7 = ^o, (ID 



where 
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is the initial value of the entropy in terms of the initial 
BH mass and Hawking temperature Tbh [4l| . 

Before proceeding, we comment briefly on the BH rate 
of absorption. The upper limit on the cross section for 
the particle absorption is [HJ 



' eft ' 



and so the accretion rate of the BH mass becomes 
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~dt 
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where 
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(16) 



is the effective BH radius for capturing particles [43| and 
e is the nearby parton energy density. The net change of 
the BH mass is therefore 



dM 
~d~t 



dM 
~df 



dM 
~dt 



(17) 
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Now, using the greybody parameters given in Table UJ it 
is easily seen that dM/dt > <^ e > 10 10 GeV/fm 3 . The 
highest earthly value of energy density of partonic matter 
will be the one created at the LHC, clhc < 500 GeV/fm 3 . 
This means that the BHs that could be produced at the 
South Pole would evaporate much too quickly to swallow 
the partons nearby. Contrary to collider experiments, 
these BHs are produced with large momentum in the lab 
system, and their decay products are swept forward with 
large Lorentz factors. 

To perform a Lorentz transformation of the evapo- 
rating BH from its rest frame S to the rest frame >S" 
of the observer (IceCube), we make use of the energy- 
momentum tensor of the outgoing Hawking radia- 
tion. In the rest frame of the BH, at a distance large 
with respect to the Schwarzschild radius of Eq. ([6|) , the 
energy-momentum tensor of the outgoing particles with 
energies uj in the range dio and with directions lying in a 
solid angle dfl is 



n^n u ■ dD 

U1 Z 4-7T 



(18) 



where iVj is defined by Eq. ([3]) as the emission rate per 
degree of freedom i of particles having energies to in the 
range dio and spin s. As noted earlier, we are working 
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with d — 4. In spherical coordinates with the BH cen- 
tered at the origin, the 4-momcntum is 



P' 



to (sin 9 cos 4>, sin 9 sin (b, cos 9, 1). 



(19) 



Integrating the (0,0) component of Eq. (fl8|) over all di- 
rections and energies, we see that T 00 is the rate at which 
the BH radiates energy per degree of freedom i of a par- 
ticle of spin s. This requirement dictated the expression 
that we wrote in Eq. flTSJ) 

Explicitly, we have for dT^ v in the BH's rest frame 
(with D = 10 and d = 4): 



dT^ = /p"^w[exp(w/T)-(-l) 2s ]- 1 



x sin < 



32tt 3 
d9 d<b du> 



(20) 



Let the black hole move at speed «bh in the z' direction 
relative to the rest frame S' of the observer, and take the 
axes of S and S' to be parallel. Then one has </>' = <fi, 



Lu(l + «BH COS ( 



(1 



^BH 



)l/2 



and 



cost/ = 



v bh + cos 9 



1 + «BH COS 9 

from which it follows that 

smO'de'dcj)' = (u;/u/) 2 sin< 
Therefore, in the rest frame of the observer 
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where 
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This can also be written using the inverse Lorentz trans- 
formation as 



T' 



(1 ~ ^1h) 1/2 

(1 - W B H COS 9') 



T, 



(26) 



where 9' is the angle between the direction of the black 
hole and of the emitted particle as measured in the ob- 
server's rest frame. Similarly, we can write 



Ld'(l — 1>BH COS#') 



(1 



BH 



)l/2 



(27) 



Putting this into the grey body factor a s (w) gives its de- 
pendence on the energy and direction of the emitted par- 
ticle as measured in the observer's rest frame. Figure [2] 
shows how particles tend to be emitted in the direction 
of motion of a moving black hole. 




FIG. 2: Distribution of Particles emitted by a BH at the origin 
and moving in the z' direction at speed 0.95c. In the BH's 
rest frame, the particle directions are uniformly distributed 
and each particle has 1 unit of energy. The length and angle 
of each line represents the energy and angle of an emitted 
particle in the observer's rest frame. Rotate the hgure about 
the z' axis for the 3-dimensional distribution. 



TABLE I: Degrees of freedom of particle species and greybody 
parameters as defined in Eqs. (0 and (|10p . 



particle's spin 


Ci 


r« 
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0.80 


0.80 
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2 


90 


0.66 


0.62 
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27 


0.60 


0.67 
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Despite the large Lorentz factors, it is easily seen from 
Eqs. (j2"Tj) , and ([2"5)) that the expected average energy 
of outgoing muons from a BH explosion is much smaller 
than the average energy of a secondary muon produced 
in a CC interaction of a neutrino with the same energy. 
Specifically, the inclusive cross section for the reaction 
fj,~~ + anything is 



dace _ 2 g| m N E v 
dy 7r 



Q 2 



x / dx{xf q (x,Q)+xf 7I (x,Q)(l-y) 2 },(28) 

where tojv is the mass of the nucleon, x = Q 2 /2miqv, 
y = v/E v , with — Q 2 the momentum transfer between 
the neutrino and muon, and v the lepton energy loss 
in the lab frame, v — E v — E^. Here, Gf = 1.16632 x 
10~ 5 GeV~ 2 is the Fermi constant, Mw is the mass of the 
W gauge boson, and f q (x, Q) and fq(x, Q) stand for com- 
binations of quark and anti-quark pdf's, respectively 4511. 
For </s > 5 TeV, the mean inelasticity is (y) < 0.25 4a ]. 
Therefore, SM events mimicking the explosion of BHs 
with So > 10 are far in the tail of the y distributuion. 
Because of this, as we show in the next section, the ra- 
tio of muon to hadronic energy deposit in a given event 
provides a clean signal to search for BHs at the South 
Pole. 



III. ICECUBE DISCOVERY REACH 

To evaluate the BH discovery reach at IceCube, one has 
to estimate the "beam luminosity," i.e. the magnitude of 
the (yet to be detected!) cosmic neutrino flux. We know 
that cosmic accelerators produce particles with energies 
in excess of 10 11 GeV (we do not know where or how [47|). 
and a neutrino beam is expected to come in association 
with these cosmic rays (48|. However, given our ignorance 
of the opacity of the sources, it is difficult to calculate the 
magnitude of the neutrino flux. The usual benchmark 
here is the so-called Waxman-Bahcall (WB) flux [49( 



~ 6.0 x lQ-^/GeV)- 2 GeV" 1 cm 



— 9 — 1 

z s sr 



(all flavours), which is derived assuming that neutrinos 
come from transparent cosmic ray sources, that the onset 
of the extragalactic component in the cosmic ray spec- 
trum is at ~ 10 10 GeV, and that there is adequate trans- 
fer of energy to pions following pp collisions. We too will 
use the WB flux to estimate the event rates necessary to 
quantify the sensitivity to BH production. To evaluate 
the sensitivity to the assumed flux, we also consider as 
an upper limit the (AARGHW) flux gS], 



3.5 x 10" 3 (i^/GeV)- 2 ' 54 GeV 



-i — ? - 
cm s 



sr 



(30) 

which is expected if extragalactic cosmic rays (from 
transparent sources) begin dominating the observed spec- 
trum at energies as low as ~ 10 86 GeV, as suggested by 



TABLE II: Expected number of BH events for M 10 = 1 TeV 
and different values of the infrared cutoff. We have taken an 
integration time of 15 yr corresponding to the lifetime of the 
experiment and used the new (old) values of F and y. The 
event rates roughly scale oc M^ 6 ^ 7 . 



^min 


Abh 


[WB] 


Abh [AARGHW] 


3 


43 


(19) 


69 (30) 


4 


34 


(15) 


43 (19) 


5 


27 


(12) 


28 (12) 


6 


22 


(9) 


20 (9) 



recent HiRes data [5 II ]. A similar enhancement in the 
neutrino flux is expected from "hidden" sources which 
are opaque to ultra-high energy cosmic rays (52j. 

IceCube is sensitive to both downward and upward 
coming cosmic neutrinos. However, to remain conser- 
vative with our statistical sample, here we select only 
downward going events. To a good approximation, the 
expected number of BH events at IceCube is given by 



A/bh = 27rn T T / dE v a{E v ) <j> v {E v ) 



(31) 



where rtx is the number of target nucleons in the effective 
volume and T is the running time of the experiment. In 
our analysis we are interested only in contained events, 
for which an accurate measurement of the inelasticty can 
be obtained. IceCube's effective volume for (background 
rejected) contained events is roug hly 1 km 3 [Hj], which 
corresponds to nx — 5.4 x 10 38 . In order to have a good 
energy resolution to determine the energy fraction in the 
muon track, we set the upper limit in Eq. (13ip well below 
detector saturation, £„ iBlax = 10 10 GeV [54(. The lower 
limit on this integral will be set so as to minimize the 
background from SM events and consequently depends 
on the infrared cutoff x m m (more on this below) . To give 
an idea of the overall picture, the total number of BHs 
expected to be produced within the lifetime of the exper- 
iment, for different "beam luminosities" and considering 
-Evmin = R) 7 GeV, are summarized in table HT1 



In the spirit of [6j , we consider the signal of BH events 
with total multiplicity N > 4 and at least one fjr 1 in the 
final state. To implement the first cut we make use of the 
average multiplicities (N) for the various particle species 
(incorporating evolution effects during Hawking radia- 
tion) summarized in the previous section. To implement 
the second cut we define the average multiplicity for any 
subset of states {s} as usual, (Nt s \) — Bi s -i(N), where 



B [s } - 



E 



6M 



Ci gi Ti 



J2i °i 9i r i 



(32) 



is the so-called "branching fraction" . Now, using the the 
parameters given in Table [fl we find (N) = 0.30 M/T 
and (AT M ±) = 0.022 (N) = 0.007 M/T. 

(N) is the average value of a Poisson distribution. If all 
species are Poisson distributed, then the sum of particles 
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in any subset is also Poisson distributed, and so N, N^± , 
and N — N,,± are all Poisson distributed, where N.,± is 
the total number of fi per event. The signal probability 
(i.e., that a given event has N^± > 1 and N > 4) is [19( 



i=0 



l 



and so the number of signal events becomes 



J dE v o(E v ) <j> v (E v ) P si 



(33) 



(34) 



The quarks and gluons emitted by the BH promptly frag- 
ment into hadrons (mostly ir^ and 7T°). For E„ > 1 TeV, 
the interaction mean free path of ir^ in ice is orders of 
magnitude smaller than the pion decay length, and so 
nearly all the hadronic energy is channeled into electro- 
magnetic modes through n° decay. The signal of such a 
hadronic/electromagnetic cascade is a bright, point like, 
source of Cerenkov light. The shower topology can be 
easily identified by the sphericity of the light pattern. 
The measurement of the radius of the lightpool mapped 
by the lattice of photomultiplier tubes determines the en- 
ergy and turns IceCube into a total absorption calorime- 
ter [55|. On the other hand, the muons emitted by the 
BH would produce a track moving outwards from the in- 
teraction vertex, providing very useful tags for the event. 

As discussed in the previous section the SM back- 
ground masking BH events are in the tail of the CC y 
distribution: for (N) > 4, the average energy of the emit- 
ted muon (after considering energy losses due to classi- 
cal radiation) is less than 20% of the incoming neutrino 
energy. Therefore, to filter the background we evaluate 
Eq. J2S]) for y > 0.8, yielding 



T y>0.8 



1.2 (E u /GeVy>- i5H pb 



(35) 



Now, substituting the CC SM cross section (with y > 0.8) 
into Eq. (|3~lj) leads to a straightforward calculation that 
shows that for E^mm = 10 s GeV the expected number 
of background events is negligible, less than 1 event in 
15 years, independently of the selected (WB/AARGHW) 
beam luminosity [5q | . The sensitivity of IceCube to probe 
the x m i n /Mio parameter space at the 3<r level (57| is sum- 
marized in Table Hill 

The expected number of background events rises with 
decreasing the low-energy cutoff. For example, for 
£v,min = 10 7 GeV and a beam luminosity at the 
AARGHW level, A/bk'd = 10. To remain conservative, 
we adopt this energy range and require a 5a excess for 
discovery. The resulting reach is shown in Fig. [3] 

The BH entropy is a measure of the validity of the semi- 
classical approximation. For x mul > 3, Sq ^> 10 yields 
small thermal fluctuations in the emission process [58| . 
Therefore, strong quantum gravity effects may be safely 



TABLE III: Sensitivity of IceCube at the 3ct level for the value 
of Mio /TeV using fiducial beam luminosities. We have taken 
an integration time of 15 yr corresponding to the lifetime of 
the experiment and used the new (old) values of F and y. 



3^ min 


Mio/TeV [WB] 


Mio/TeV [AARGHW] 


3 


1.5 


(1.2) 


1.5 


(1.2) 


5 


1.3 


(1.1) 


1.3 


(1.1) 


7 


1.2 


(1.0) 


1.2 


(1.0) 


9 


1.1 


(1.0) 


1.1 


(0.9) 



> 



4.5 



3.5 



2.5 



0.5 



LH( 



3.5 



4.5 



FIG. 3: IceCube discovery reach of quantum BHs, assuming 
the new estimates of F and y given in For comparison 

the LHC discovery reach, assuming a cumulative integrated 
luminosity of 1 ab _1 over the life of the collider, is also shown. 
The LHC discovery reach has been obtained by scaling up the 
results given in 19], to account for the new values of y and 
F. 



neglected in this "energy regime," Mj^/MiQ > 3. More- 
over, gravitational effects due to brane back- reaction 
are expected to be insignificant for Mbh well beyond 
the brane tension, which is presumably on the order of 
Mio [lil • As noted in the previous section, string theory 
provides a more complete picture for Mbh close to M w . 
In string theory, the ultimate fate of the BH is determined 
by the string ^ BH correspondance principle: when the 
Schwarzschild radius of the BH shrinks to the fundamen- 
tal string length l s ^> 10~ 19 m an adiabatic transition oc- 
curs to a massive superstring mode |32| . Subsequent en- 
ergy loss continues as thermal radiation at the unchang- 
ing Hagedorn temperature 59]. The continuity of the 
cross section at the correspondance point, parametrically 
in both the energy and the string coupling, provides inde- 
pendent support for this picture [201 ] . In the perturbative 
string regime, however, the parton-parton cross section 
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contains the Chan Paton factors which control the pro- 
jection of the initial state onto the string spectrum. In 
general, this projection is not uniquely determined by 
the low-lying particle spectrum, yielding one or more ar- 
bitrary constants [601 ]. Interestingly, the parton-parton 
cross section derived in [2(| from the Virasoro-Shapiro 
amplitude leads to an enhancement of the predicted BH 
cross section for 1 < x m j n < 3 [61(. This makes it plau- 
sible to adopt the BH cross section as a lower bound of 
the real Planckian cross section. However, it is important 
to stress that the proposed signal to search for BHs at 
IceCube strongly depends on the large density of quan- 
tum mechanical states of BHs with So ^S> 10, and conse- 
quently it is only valid in the semi-classical regime. The 
LHC signal (hard photons and leptons in the collision 
rest frame), however, can be mantained under plausible 
hypotheses on the superstring decay modes. With this in 
mind, the LHC discovery reach for the quantum regime 
shown in Fig. [3] may be taken as a lower bound, derived 
on the assumption that the BH cross section provides a 
lower limit on the string cross section. 

IV. CONCLUSIONS 

In this work we have reviewed the possibility of search- 
ing for BHs using cosmic neutrino interactions in the 
Antarctic ice. We have shown that the ability of Ice- 
Cube to accurately measure the inelasticity distribution 
of events provides a unique discriminator for SM back- 
ground rejection [§2], allowing extremely sensitive probes 
of TeV-scale BH production. In the optimistic case that 
the neutrino flux is at the WB level, IceCube has a sub- 
stantial discovery potential for BHs, well in the semi- 
classical regime where M^ 1 > 3Afio- The statistically 
significant 5tr excess extends up to M 10 ~ 1.3 TeV. 

A point worth noting at this juncture: In assessing 
the discovery potential of IceCube we have performed 
a Lorentz transformation of the evaporating BH from 
its rest frame to the lab frame, so as to obtain the 
boosted energy spectrum of the outgoing Hawking ra- 
diation. Such a spectrum is also relevant for the study 
of near-extremal supersymmetric BHs, which could be 
copiously produced at the LHC. These BHs can be as- 
sociated with the supermassive string modes expected to 
populate the quantum regime [63| ■ Note that in order to 
mantain the configuration of the initial state with "zero 
supersymmetry," the central charge conservation would 
force these BHs to be pair produced with non-zero trans- 
verse momentum (pt), travelling along the beam pipe. 
Therefore, the prompt decay of the BHs would produce 
a startlingly clean signal that should have very few back- 
grounds: their decay products (that of course may in- 
clude sparticles [H)]), would trigger high multiplicity fire- 
works with boosted spherical shape collimated into back- 
to-back pencil beams. A crude estimate of the event rates 
can be obtained from the analysis of u ij — > BH + others" 



subprocesses (65[. The results are encouraging: the pro- 
duction rate of BHs with large pt > 500 GeV (which 
is balanced by the momentum of an energetic parton) 
would still be large enough for detection at LHC (assum- 
ing an integrated luminosity 100 fb _1 ). Of course the 
energy requirement for BH pair production would yield 
an additional suppression factor and its proper consider- 
ation requires a Monte Carlo simulation. Therefore, we 
strongly recommend to include in future versions of the 
BH event generators CHARYBDIS and catfish [67| 
a detailed treatment of production and evaporation of 
near-extremal supersymmetric BHs. 

In closing, it is important to stress that this analy- 
sis is meant to investigate the underlying principles and 
does not account for the details of the detector response. 
Though the devil is generally in the details, we believe 
our estimate of the IceCube discovery reach is conserva- 
tive, as we have not included the r-channel which has 
the potential of nearly doubling our signal event back- 
ground. Specifically, r leptons emitted by the BH may 
decay in flight inside the instrumented volume after es- 
caping the electromagnetic shower bubble, thereby trig- 
gering a "second bang [68| ." The ratio of the first to sec- 
ond bang fractional energy provides a clean direct signal 
of BH production, and like the muon channel is "inde- 
pendent" of the absolute neutrino flux. The inclusion of 
the tau channel, however, requires a full blown Monte 
Carlo simulation to properly determine the acceptance 
for such events, in which the triggering probability for 
the second bang depends on Tbh and the corresponding 
tau decay length. 

In summary, over the next few years high-statistics 
high-energy precision data will be collected at the LHC. 
In addition, the IceCube neutrino telescope is coming on 
line with complementary information at ultra-high ener- 
gies. This new arsenal of data will certainly provide an 
ideal testing ground for TeV-scale BH production, and, 
at the same time, a unique opportunity to view similar 
physics from two different points of view. Should the 
LHC find evidence of BHs a bit outside the range acces- 
sible for the baseline IceCube design, the ideas discussed 
in this paper could constitute another compelling reason 
for pursuing HyperCube [69| . 
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